A model for giant flares in soft gamma repeaters 
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We argue that giant flares in SGRs can be associated to the core conversion of an isolated neutron 
star having a subcritical magnetic field ~ 10 12 G and a fallback disk around it. We show that, in a 
timescale of < 10 5 yrs, accretion from the fallback disk can increase the mass of the central object 
up to the critical mass for the conversion of the core of the star into quark matter. A small fraction 
of the neutrino-antineutrino emission from the just-converted quark- matter hot core annihilates into 
e + e~ pairs above the neutron star surface originating the gamma emission of the spike while the 
further cooling of the heated neutron star envelope originates the tail of the burst. We show that 
several characteristics of the giant flare of the SGR 1806-20 of 27 December 2004 (spike and tail 
energies, timescales, and spectra) can be explained by this mechanism. 



Soft 7-ray repeaters (SGRs) are persistent X-ray emit- 
ters that sporadically emit short bursts of soft 7-rays. 
In the quiescent state they have an X-ray luminosity of 
~ 10 35 erg/s, while during the short 7-bursts they re- 
lease up to 10 42 erg/s in episodes of about 0.1 s. Only 
four SGRs are known at present, 3 in our Galaxy, and 
one in the Large Magellanic Cloud. Exceptionally, three 
of them have emitted very energetic giant flares which 
commenced with brief 7-ray spikes of ~ 0.2 s, followed 
by tails lasting hundreds of seconds. Hard spectra (up 
to 1 MeV) were observed during the spike and the hard 
X-ray emission of the tail gradually faded modulated at 
the neutron star (NS) rotation period. 

The powerful giant flare of the SGR 1806-20 of De- 
cember 27, 2004 is particularly important because it was 
observed by many satellites [lj, a a a a. The duration 
of the initial spike was ~ 0.25 s. For the first 500 ms 
the emission was so strong as to saturate almost all the 
detectors on 7-rays satellites, though there are unique 
measurements of the initial flare activity The lumi- 
nosity of the spike was ~ 10 47 erg/s (assuming a distance 
of 15 kpc and isotropic radiation), hundreds of times 
brighter than the giant flares previously observed from 
other SGRs. A long tail lasting ~ 400 s followed the 
initial spike with a thermal component with decreasing 
temperatures in the ~ few keV range [l| . There is a clear 
modulation with the same period of 7.56s already known 
from previous studies of the quiescent emission. 

While several characteristics of SGRs are often ex- 
plained in terms of the magnetar model, assuming that 
the object is a NS with an unusually strong magnetic 
field (B ~ 10 15 G) 0, there is alternative theoretical 
work trying to explain the behavior of SGRs and Anoma- 
lous X-ray Pulsars (AXPs) by assuming the existence of 
fallback disks around isolated NSs with subcritical mag- 
netic fields of ~ 10 12 G (@,H 0,0 and Refs. therein). 
NSs are believed to be born from core— collapse super- 



novae with masses near the Chandrasekhar limit of the 
Fe core (~ 1.4M©). During the supernova explosion, 
masses typically as large as ~ 0.2M Q may fall back [ll[ 
and although most of this material will be directly ac- 
creted onto the NS, part of it can form an accretion disk 
within a few hours of the initial explosion 0. In fact, a 
fallback disk around an AXP has been recently detected 
for the first time [p| ■ In this work we will show that ac- 
cretion from a fallback disk can increase the mass of the 
NS up to the critical value for the conversion of the core 
into quark matter and this will trigger a giant 7-fl.are. 

The disk of SGRs: The inner disk radius R m of 
a fallback disk can be approximately evaluated from 
the balance between the stellar magnetosphere stress 
and the ram pressure of the inflow matter R m ~ 
[y 2 /{2GM) 1 / 2 M] 2 / 1 , where M is the mass of the NS, 
/i = R 3 B its magnetic moment, and M is the mass in- 
flow rate through the disk. If the NS magnetosphere ro- 
tates faster than the inner disk, i.e., if R m > R c , where 
R c = [GM/fl 2 ] 1 / 3 is the corotation radius (defined as the 
radius where the disk rotates with the angular velocity of 
the stellar magnetosphere), then the system is in the so 
called propeller regime in which there is an angular mo- 
mentum flux from the NS to the disk (that will lead to 
strong stellar spin-down as required by the observations 
of SGRs and AXPs). In the case that the NS rotates only 
a little faster than the inner disk (R m J> R c ), a signifi- 
cant fraction of the mass lost by the disk can accrete onto 
the NS, while at the same time the NS remains spinning 
down in a " tracking" regime @, ■ In recent work, Eksi 
et al. Q identified these different regimes in a stellar 
period P versus mass inflow rate M diagram. 

As in a S > we shall assume here that SGRs are NSs 
with moderate .B-fields that spin down to P » 5 — 8 s, 
as observed, in timescales <; 10 5 yrs due to a fallback 
disk. Notice that, if attributed to electron resonance, 
the observed 5 keV cyclotron lines of 1806-20 12j suggest 
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B ~ 5 x 10 12 G. After a brief accretion phase that follows 
its formation, the disk enters the propeller (spindown) 
regime and later, it reaches the tracking regime. During 
this phase, the system will lie at the boundary between 
the propeller and the X-rays luminous accretor (shaded) 
zone of the P — M diagram of Ref . p| . According to that 
diagram, a NS with B ~ 5 x 10 12 G can enter the tracking 
phase only if its period is J> 0.2 s and M < 10 20 g/s. Once 
reaching this point, we assume that the disk will evolve 
along the tracking path until the present observed period 
of 7.56 s (for SGR 1806-20) which corresponds to a disk 
mass inflow rate M ~ 5 x 10 16 g/s @|. 

Quasi-steady luminosity and frequent soft "/-flares: 
The disk mass accretion rate onto the NS during the 
tracking regime may consist of two contributions. One 
of them is a quasi-steadily evolving sub-Eddington com- 
ponent, M qs , whose value is of the order of M above. 
In the the framework of the fallback disk model, the 
persistent X-ray luminosity of ~ 10 35 erg/s produced 
during the more quiescent states of the SGR 1806-20 
is attributed to this component. The maximum lumi- 
nosity that a NS being presently powered by an ac- 
cretion rate M qB ~ M ~ 5 x 10 16 .g/s can produce is 
L ~ GMM/R ~ 2.7 x 10 36 erg/s, so that 4% of the 
accretion power will be sufficient to explain the observed 
persistent X-ray luminosity. The other disk accretion 
rate component may be associated to the observed soft 
7-ray pulses of SGRs with luminosities L ~ 10 42 erg/s. 
A number of mechanisms have already been suggested 
to explain these pulses. Among them, a hybrid model 
(proposed, e.g., in [9() assumes the existence of local 
very strong multipole components of the -B-field with 
B ~ 10 14 - 10 15 G that act on the NS crust to pro- 
duce the flares in a similar way to the magnetar models, 
while the large scale dipole component is ~ 10 12 — 10 
G, as required in the fallback disk scenario. Alterna- 
tively, we could speculate here that a sporadic super- 
Eddington disk mass accretion rate M sp ~ 5 x 10 21 g/s 
onto the NS surface could produce the observed flares 
with L ~ GMM/R ~ 10 42 erg/s (cf. Q). Violent mag- 
netic reconnection events occurring between the NS mag- 
netosphere and the inner disk lines could t rigg er these 
events with the observed rise times ~ 10 ms [15| . 

Accretion and conversion of the core: The total ac- 
creted mass onto the NS during the source lifetime will 
be the sum of both accreting components above, and may 
be written as J M qs dt+ Nr sp x M sp , where N is the num- 
ber of soft bursts during the object's lifetime (several 
hundreds [6j) and the integral is taken along the tracking 
path between M qs ~ 10 20 g/s and 5 x 10 16 g/s. Using a 
self-similar evolution for M qs oc t~ a , with a ~ 1.2 Q, we 
obtain a total accreted mass AM ~ 1O _3 M in a time 
interval <J 10 5 yr, most of which comes from the first 
component. This increase of the NS mass is essential 
for triggering the conversion of the NS core into quark 



matter (QM). 

In this paper we are not assuming the QM to be ab- 
solutely stable (i.e. we assume that it has an energy per 
baryon at zero pressure and temperature that is larger 
than the neutron mass). Thus, stars containing QM are 
hybrid, i.e. contain /3-stable QM only at their interior 
and not up to the surface (as in strange stars). The crit- 
ical mass M cr for the formation of a QM core in a cold 
hadronic NS has been calculated recently exploring the 
effect of surface tension and color superconductivity [lj| . 
If the parameters of the hadronic and QM equations of 
state are such that QM is not absolutely stable, M cr is 
very close to (but smaller than) the maximum mass of 
hadronic stars for a large range of the parameter space 
[l6j . We emphasize that the conclusions of the present 
paper are not sensitive to the exact value of M cr provided 
that it is larger than the NS mass short after the super- 
nova explosion and smaller than the maximum mass of 
NSs. Thus, if right after the initial supernova fallback, 
the NS acquires a mass that is not too far from M cr , then 
the small accreted AM above will be sufficient to trigger 
the conversion of the core of the NS into QM at anytime 
after the beginning of the tracking phase. 

Energy release from the conversion: The conversion 
into QM of a large part of the core releases a large amount 
of energy AE that can be estimated by the difference be- 
tween the gravitational mass of the initial and final con- 
figurations (both having the same baryonic mass) [lj|. 
As shown recently [lf|, AE is ~ 10 53 - 10 53 6 ergs for the 
conversion of a hadronic star having the critical mass. 
We assume that most of this energy is thermalized inside 
the quark matter core, in the form of a trapped gas of 
neutrinos and antineutrinos of all flavors iviVi pairs). To 
be sure, a small fraction can be used to excite vibrating 
modes 19], but we shall neglect this mechanical energy 
in the rest of the paper. Then, the initial temperature of 
the core is estimated from AE ~ C q T, where C q is the 
heat capacity of the quark core, resulting T\\ ~ 1 [23| . 

Neutrino annihilation: Once they leave the opaque and 
hot inner core, the ViVi pairs can annihilate into e + e~ 
pairs in the outer layers of the NS, and also above the 
NS surface. Annihilation within the NS will heat mat- 
ter significantly while above the NS surface it will cre- 
ate a 7-burst through e + e~ — > 7. For simplicity, we 
shall assume a blackbody spectrum for the vv's in the 
core. Also, for annihilation above the NS surface, we 
can neglect blocking effects in the phase spaces of 
and e + . Thus, the "unblocked" local energy deposition 
rate at a radial position r due to z/jOj — > e + e _ can be 
written as Q™+ b e - = AT^xlx], being T v \\ the temper- 
ature at the ^-sphere (i.e. the last scattering surface of 
the vis's), A = 1.28 x 10 34 erg cm _3 s _1 for i/ e i/ e 's, and 
A = 2.7 x 10 33 erg cm _3 s _1 for v^D^s and z/ T £v's [l7j ]. 
The angular factor x[x] — g(l — x) 4 (5 + 4x + x 2 ) depends 
on x(r) = (1 — R^/r 2 ) 1 / 2 , being R v the J/-sphere radius. 
X[x] is responsible for the rapid decrease in Q^+ b e - for 
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increasing r, as the interacting neutrinos become more 
collinear. The total luminosity injected above the NS 
surface is determined by L e+e - = f R ' Q^+ h e _Aitr 2 dr = 

£,AT^R Q where £ = J^ (R) 4ttx{1 - x 2 )- b / 2 x[x]dx is a 
function of the ratio between R v (which we assume to be 
coincident with the quark-matter core radius Rq) and 
the NS radius R. Within the NS, electrons are degen- 
erate and there is a strong blocking effect over Q e + e -. 
This can be accounted for by applying an average block- 
ing factor to the unblocked results [l7| • Then, one writes 
Q e+e - = BeQ^ b e _, where B e ~ l/[l + exp(( Me -T v )/T)], 
fi e is the e~ chemical potential at the annihilation point, 
and T the corresponding temperature jl7j |. 

Cooling evolution of the NS: The conversion of the in- 
ner core region into quark matter is expected to be very 
fast (~ 10" 5 s 18]). We assume for simplicity that af- 
ter the conversion, the NS settles reasonably fast into a 
nearly hydrostatic configuration. Though this assump- 
tion is not strictly correct, particularly during the first 
seconds of the evolution, we shall see below that the char- 
acteristics of the 7-ray tail over 400 s are determined by 
a thin shell near the surface of the NS which is much 
less sensitive to dynamical rearrangement than the core 
region. As we will also see, it is the neutrinos from the 
core that produce the spike. Even if the core is subject 
to dynamical variations these will hardly affect the order 
of magnitude of the released energy and the ^-emission 
timescale in comparison to the hydrostatic case (cf. 19]). 
Therefore, soon after the conversion we adopt a simpli- 
fied model in which the NS (with M — M cr ) has a hot 
quark matter inner-core with Rq m 0.6R in equilibrium 
with a trapped vv gas at T v \\ « 1 surrounded by an 
initially cold outer-core (composed of free nucleons, elec- 
trons and muons) and a crust (composed of nuclei, free 
neutrons and electrons). We assume there is no signif- 
icant lepton number gradient, since most of the z/s are 
produced in pairs, so that the lepton-number-transport 
equations will not dominate the evolution. Thus the cool- 
ing of the quark core is described by the energy-transport 
equation: 
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where the neutrino energy density is e„ = 3 x |aT 4 , a = 
7r 2 /c|j/15c 3 ?i 3 , and c q is the specific heat of quark matter. 
The neutrino mean free path A « 1.3 x 10 3 cm/T 1 3 1 in the 
hot inner-core is very small and the vv's are released in a 
diffusion timescale, due to scattering with free quarks. At 
the surface of the quark core we assume a free streaming 
regime into vacuum. Solving numerically this equation 
we find the time evolution of the temperature T v at the 
^-sphere. T v is then used to evaluate the luminosities 
L e + e - given above and L v = 47ri?g x 3 \osbT^ (Fig. 1). 

The cold outer-core and the crust are essentially trans- 
parent to vv's and therefore they will hardly heat due 
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FIG. 1: Left: The emission spectrum according to the present 
model (dashed line) compared with the observed spectrum of 
SGR 1806-20 (dots with errors included) given in Right: 
Luminosity of the core neutrino emission and of the e + e~ 
injection above the NS surface (~ L 7 of the spike) during the 
first few seconds after the conversion. 



to the emission of the hot inner-core. Also, since e~'s 
are strongly degenerate, heating of a given layer through 
vv — > e + e" will be significant only if the mean energy 
of the vv's (which is of the order of that of the outgo- 
ing e + e") is larger than the Fermi energy of the e~'s 
at that layer, i.e. ksT v ;> fi e . Since T v <; 10 MeV, 
this is fulfilled only in the outer-crust /envelope. Using a 
simple relation between the depth z measured downward 
from the NS surface and fi e [20j, z w 25m(/i e /r7i e — 1) s» 
25m[X'„/m e — 1], we find that initially only the outermost 
~ 500 meters will be significantly heated by this process. 
Since the calculations show that T v (t) decreases expo- 
nentially (Fig. 1), heating ceases in less than ~ 1 s and 
thereafter, the envelope evolves decoupled from the core. 
The cooling of the envelope can then be described in a 
plane parallel approximation by: 
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where the specific heat is approximately that of a Dulong- 
Petit solid and a degenerate relativistic electron gas, 
c v w (4.5xl0 16 pio + 3.9xl0 16 pio 3 r 9 )ergcm- 3 K" 1 . The 
thermal conductivity is K = lQassT 3 / (3/?k) with kT 1 — 
K cond + K rad where the conductive opacity IS Kcond — 
3.6 x 10" 2 pg _1 ' 2 T 9 1 - 4 cm 2 g _1 and the radiative opacity is 
Krad = 9.6 x 10" 2 pg 1 r 9 ~ 1 - 2 cm 2 g~ 1 (within a factor of 
~ 3) 2l|. The density profile is given by p6 ~ (z/10 4 cm) 3 
[20I ] and the effective surface is taken where the opti- 
cal depth is ~ 2/3. For simplicity, the ^-cooling sink is 
given approximately by Q v = 1.03 x 10 15 To 6 erg cm" 3 s" 1 



for Tg > 1.3 and 10~ 2 £ p 6 ^ 10~ 2 |21(. The energy 
per unit volume injected by core-z/s in the envelope is 
U(z) ~ j B e ATl lx {t)x[x]dt. Assuming that this energy 
is thermalized the initial temperature profile is given by 
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~ U (z) / c v (z) > 10 10 K, but due to the efficient i^-cooling 
sink it falls to values < 10 10 K in less than 1 s. Thus the 
later evolution of the envelope is rather insensitive to the 
initial energy injection. 

Results: In the right panel of Fig. 1, we show the 
time evolution of the core neutrino luminosity and the 
luminosity L e + e - of the pairs that annihilate above the 
NS's surface. The timescale at which L e + e - decays to 
~ 1% of its initial value is ~ 0.2 s, just as observed for 
the spike of giant flares in SGRs (this is a robust result 
that naturally arises also in models of GRBs involving 
the burning of quark stars; e.g. [22j). We also note that 
the efficiency of the vv — > e + e~ conversion is strongly 
temperature dependent, so that with variations within a 
factor of ^ 2 — 3 in the central temperature of the object, 
a wide range of observed spike-luminosities (10 44 — 10 47 
erg/s) can be explained. The spectrum of the spike has 
been calculated assuming that all pairs are converted into 
7-rays (i.e. L e+e - fa L 7 ). This results T 7 (f) oc T*(*(t) 
which reproduces the observed best-fit blackbody tem- 
perature T splke = 175 ± 25 keV of SGR 1806-20 Q for an 
initial core temperature T v \\ ps 0.5 (left panel of Fig. 1). 
The luminosity L e + e - fades very fast and then is over- 
whelmed by photon radiation from the hot NS envelope. 
The surface temperature evolves according to Eq. (2) 
producing the tail emission as shown in Fig. 1. 

Discussion: The model just described can explain the 
energy scale, spectrum and timescale of both the spike 
and the tail of giant flares in SGRs. Several features 
can be used to observationally distinguish between this 
and the magnetar model. For example, the concomi- 
tant detection of low energy neutrinos and anti-neutrinos 
(E u ps 10 MeV) with the 7-emission will give strong sup- 
port to our model. Unlike in the magnetar model, a 
huge amount of energy stored in the core is converted in 
7-emission with just a tiny efficiency The eventual ob- 
servation of a giant flare with energy much larger than 
10 47 erg will also give support to this model, while it 
would be hard to explain within the magnetar scenario. 
Also, all four SGRs are possibly associated with super- 
nova remnants (SGR 1900+14 / G42. 8+0.6, SGR 0526- 
66 / N49, SGR 1806-20 / G10.0-0.3, and SGR 1627-41 / 
G337. 0-0.1) which makes it plausible that they still have 
a fallback disk. The recent identification of a fallback 
disk around the AXP 4U 0142+61 [l]| points to this di- 
rection, provided that SGRs and AXPs belong to the 
same population. The release of energy in the core and 
the NS-structure rearrangement after phase transition, 
as well as the sporadic impact of accreted matter, may 
cause surface stresses and quakes just as those that in 
the magnetar scenario are caused by magnetic stresses. 
Therefore, quasi-periodic-oscillations (QPOs), for exam- 
ple, could be likewise explained by crust fracturing and 
release of elastic energy. Finally, since AXPs seem to 



have bursts but not giant flares, we can interpret them 
as accreting NSs whose mass is far from M cr . In this 
sense, the number of AXPs should be larger than that 
of SGRs. Also, our model predicts that all SGRs that 
had a giant flare should have essentially the same mass, 
i.e. M cr . Techniques for determination of NS- masses 
are rapidly improving and eventual measurements could 
probe the value of M cr . If so, this could also provide key 
information for the equation of state of dense matter. 
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FAPESP and CNPq. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 

[8] 

[9] 
[10] 

[11] 

[12] 
[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 



[20] 
[21] 

[22] 

[23] 



Electronic address: german.lugones@ufabc.edu.br 
K. Hurley et al. Nature 434, 1098 (2005) - 
Schwartz S. J. et al. Astrophys. J. 627, L129 (2005) 
S. Mereghetti et al. Astrophys. J. 624, L105 (2005) 
T. Terasawa et al. Nature 434, 1110 (2005). 

Palmer, D. M., Nature 434, 1107 (2005). 

P. M. Woods & C. Thompson, arXiv: |astro-ph/0406133 
P. Chatterjee, L. Hernquist, and R. Narayan; Astrophys. 
J. 534, 373, (2000) 

K. Y. Eksi, L. Hernquist, and R. Narayan; Astrophys. J. 
623, L41, (2005) 

K. Y. Eksi and M. A. Alpar, Astrophys. J. 599, 450 (2003) 
U. Ertan et al.; Astrophys. J. 657, 441 (2007) 
C. L. Fryer, S. A. Colgate and P. A. Pinto; Astrophys. J. 
511, 885 (1999) 

A. I. Ibrahim et al., Astroph. & Sp. Sci. 308, 43 (2007) 
Z. Wang, D. Chakrabarty and D. L. Kaplan, Nature 440, 
772 (2006) 

H. Mosquera Cuesta et al, Phys. Rev. Lett. 80, 2988 
(1998) 

E.M. de Gouveia Dal Pino, A. Lazarian, Astron. & As- 
trophys., 441, 845 (2005) 

I. Bombaci, G. Lugones and I. Vidana, Astron. & Astro- 
phys., 462, 1017 (2007) 

J. Cooperstein, L.J. van den Horn & E. Baron, Astro- 
phys. J. 321, L129 (1987); Astrophys. J. 309, 653 (1986). 
G. Lugones, O. G. Benvenuto & H. Vucetich, Phys. Rev. 
D50, 6100 (1994) 

G. Marranghello, C.Z. Vasconcellos and J. A. de Frcitas 
Pacheco Phys. Rev. D 66, 064027 (2002), G. Miniutti, 
J. A. Pons, E. Berti, L. Gualtieri and V. Ferrari, MNRAS 
338, 389 (2003). 

L. Hernquist & J.H. Applegate, Astrophys. J. 287, 244 
(1984) 

E.H. Gudmundsson et al. Astrophys. J. 272, 286 (1983); 
Y. Lyubarsky et al. Astrophys. J. 580, L69 (2002); D. 
Eichler & A.F. Cheng, Astrophys. J. 336, 360 (1989). 
P. Haensel, B. Paczynski and P. Amsterdamski, Astro- 
phys. J. 375, 209 (1991); G. Lugones et al., Astrophys. J. 
581, L101 (2002) 

We shall use the notation Y n = Y/10 n [Y], where n is an 
integer, and [Y] stands for the units of the quantity Y. 
We use cgs units and K for the temperature. 



